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Abstract 

^ ■ We confront the latest HI and ZEUS data on diffractive dijet photoproduction 

with next-to-leading order QCD predictions in order to determine whether a 
^ \ rapidity gap survival probability of less than one is supported by the data. We 

<^ ■ find evidence for this hypothesis when assuming global factorization break- 

^ , ing for both the direct and resolved photon contributions, in which case the 

fN| \ survival probability would have to be £y 6 *-dependent, and for the resolved or 

in addition the related direct initial-state singular contribution only, where it 
^ . would be independent of E^. 

1 Introduction 

& . 

The central problem for hard diffractive scattering processes, characterized by a large rapidity gap in 
high-energy collisions, is whether they can be factorized into non-perturbative diffractive parton density 
functions (PDFs) of a colorless object (e.g. a pomeron) and perturbatively calculable partonic cross 
sections. This concept is believed to hold for the scattering of point-like electromagnetic probes off 
a hadronic target, such as deep-inelastic scattering (DIS) or direct photoproduction [1], but has been 
shown to fail for purely hadronic collisions [1,2]. Factorization is thus expected to fail also in resolved 
photoproduction, where the photon first dissolves into partonic constituents, before these scatter off the 
hadronic target. The separation of the two types of photoproduction processes is, however, a leading 
order (LO) concept. At next-to-leading order (NLO) of perturbative QCD, they are closely connected 
by an initial-state (IS) singularity originating from the splitting 7 — > qq (for a review see [3]), which 
. may play a role in the way factorization breaks down in diffractive photoproduction [4]. The breaking 

of the resolved photoproduction component only leads to a dependence of the predicted cross section 
on the factorization scale M 7 [4]. Since this M 1 -dependence is unphysical, it must be remedied also 
for the factorization breaking of the resolved part of the cross section, e.g. by modifying the IS singular 
direct part. A proposal how to achieve this has been worked out in [4] and has been reviewed already 
in the proceedings of the workshop on HERA and the LHC of 2004-2005 [5] (see also [6]). Since from 
a theoretical point of view only the suppression of the resolved or in addition the IS singular direct 
component [4] is viable, it is an interesting question whether the diffractive dijet photoproduction data 
show breaking of the factorization, how large the suppression in comparison to no breaking will be, and 
whether the breaking occurs in all components or just in the resolved plus direct IS component. The 
value of the suppression factor or survival probability can then be compared to theoretical predictions [7] 
and to the survival probability observed in jet production in pp collisions at the Tevatron [2] and will be 
of interest for similar diffractive processes at the LHC. 

Since 2005 no further developments occurred on the theoretical side. On the experimental side, 
however, the final diffractive PDFs (DPDFs), which have been determined from the inclusive measure- 
ments of the diffractive structure function by the HI collaboration, have been published [8]. Also 
both collaborations at HERA, HI and ZEUS, have now published their final experimental data of the 
cross sections for diffractive dijet photoproduction [9, 10]. Whereas HI confirm in [9] their earlier 
findings based on the analysis of preliminary data and preliminary DPDFs, the authors of the ZEUS 
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analysis [10] reached somewhat different conclusions from their analysis. Specifically, the HI collab- 
oration [9] obtained a global suppression of their measured cross sections as compared to the NLO 
calculations. In this comparison [9], the survival probability is R = 0.5, independent of the DPDFs fit 
used, i.e. fit A or B in Ref. [8]. In addition they concluded that the assumption that the direct cross section 
obeys factorization is strongly disfavored by their analysis. The ZEUS collaboration, on the other hand, 
concluded from their analysis [10], that, within the large uncertainties of the NLO calculations, their data 
are compatible with the QCD calculations, i.e. that no suppression would be present. 

Due to these somewhat inconsistent results we made a new effort [11] to analyze the HI [9] and 
the ZEUS [10] data, following more or less the same strategy as in our earlier work [12, 13] on the 
basis of the NLO program of [12, 13] and the new DPDFs sets of Ref. [8]. The HI and the ZEUS 
dijet data cannot be compared directly, since they have different kinematic cuts. In particular, in the 
HI measurements [9] E 3 ^ tl{2) > 5 (4) GeV and x P < 0.03, and in the ZEUS measurements [10] 
E^ 2 ) > 7.5 (6.5) GeV and xjp < 0.025 (these and all other variables used in this review are defined 
in [11-13] and in the corresponding experimental contribution in these proceedings). It is clear that in 
order to establish a global suppression, i.e. an equal suppression of the direct and the resolved cross 
section, the absolute normalization and not so much the shape of the measured cross section is very 
important. This normalization depends on the applied kinematic cuts. Of course, the same cuts must 
be applied to the NLO cross section calculation. In case of a resolved suppression only, the suppression 
depends on the normalization of the cross sections, but also on the shape of some (in particular the i*, 
E 3 ^ a , M12, and fji ets ) distributions, and will automatically be smaller at large E^ n [3]. Distributions 
in xjp and y (or W) are not sensitive to the suppression mechanism. The distribution in zp, on the other 
hand, is sensitive to the functional behavior of the DPDFs, in particular of the gluon at large fractional 
momenta. 

In the meantime, the HI collaboration made an effort to put more light into the somewhat contra- 
dictory conclusions of the HI [9] and ZEUS [10] collaborations by performing a new analysis of their 
data, now with increased luminosity, with the same kinematic cuts as in [9], i.e. the low-E^ 6 * 1 cut, and 
the high-E^T 1 cut as in the ZEUS analysis [10]. The results have been presented at DIS 2008 [14] and 
will be published soon. We have performed a new study of these HI [14] and ZEUS data [10] to show 
more clearly the differences between the three data sets [15]. In this contribution we shall show a selec- 
tion of these comparisons. The emphasis in these comparisons will be, how large the survival probability 
of the diffractive dijet cross section will be globally and whether the model with resolved suppression 
only will also describe the data in a satisfactory way. In section 2 we show the comparison with the HI 
data [14] and in section 3 with the ZEUS data [10]. Section 4 contains our conclusions. 



2 Comparison with recent HI data 

The recent HI data for diffractive photoproduction of dijets [14] have several advantages as compared 
to the earlier HI [9] and ZEUS [10] analyses. First, the integrated luminosity is three times higher 
than in the previous HI analysis [9] comparable to the luminosity in the ZEUS analysis [10]. Second, 
HI took data with low-E^ 6 * and high-E 3 ^ 1 cuts, which allows for a comparison of [9] with [10]. The 
exact two kinematic ranges are given in [14]. The ranges for the low-^ e * cuts are as in the previous 
HI analysis [9] and for the high-E^ 6 ' cuts are chosen as in the ZEUS analysis with two exceptions. In 
the ZEUS analysis the maximal cut on Q 2 is larger and the data are taken in an extended y-range. The 
definition of the various variables can be found in the HI and ZEUS publications [9, 10]. Very important 
is the cut on xjp. It is kept small in both analyses in order for the pomeron exchange to be dominant. 
In the experimental analysis as well as in the NLO calculations, jets are defined with the inclusive hp- 
cluster algorithm [16, 17] in the laboratory frame. At least two jets are required with the respective cuts 
on E^ n and E 3 ^ t2 , where E 3 ^ 11 ^ refers to the jet with the largest (second largest) £^f*. 

Before we confront the calculated cross sections with the experimental data, we correct them for 
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Fig. 1: Differential cross sections for diffractive dijet photoproduction as measured by HI with \ow-Ei^ cuts and compared to 
NLO QCD without (R = 1) and with (R — 0.46) global suppression (color online). 
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Fig. 2: Differential cross sections for diffractive dijet photoproduction as measured by HI with \ow-E^ cuts and compared to 
NLO QCD with global, resolved, and resolved/direct-IS suppression. 



hadronization effects. The hadronization corrections are calculated by means of the LO RAPGAP Monte 
Carlo generator. The factors for the transformation of jets made up of stable hadrons to parton jets were 
supplied by the HI collaboration [14]. Our calculations are done with the 'HI 2006 fit B' [8] DPDFs, 
since they give smaller diffractive dijet cross sections than with the 'HI 2006 fit A. We then take nf = 4 
with = 0.347 GeV, which corresponds to the value used in the DPDFs 'HI 2006 fit A, B' [8]. For 
the photon PDFs we have chosen the NLO GRV parameterization transformed to the MS scheme [18]. 

As it is clear from the discussion of the various preliminary analyses of the HI and ZEUS collab- 
orations, there are two questions which we would like to answer from the comparison with the recent HI 
and the ZEUS data. The first question is whether a suppression factor, which differs substantially from 
one, is needed to describe the data. The second question is whether the data are also consistent with a 
suppression factor applied to the resolved cross section only. For both suppression models it is also of 
interest whether the resulting suppression factors are universal, i.e. whether they are independent of the 
kinematic variables of the process. To give an answer to these two questions we calculated first the cross 
sections with no suppression factor (R = 1 in the following figures) with a theoretical error obtained 
from varying the common scale of renormalization and factorization by factors of 0.5 and 2 around the 
default value (highest E^ 1 ). In a second step we show the results for the same differential cross sections 
with a global suppression factor, adjusted to da/dE^ 11 at the smallest E^ a -bm. As in the experimental 
analysis [14], we consider the differential cross sections in the variables x° bs , Zp s , log 10 (xjp), E^ tl , 
Mx, M\2, rf ets , \^.rf ets \ and W [15]. Here we show only a selection, i.e. the cross sections as a func- 



ep -> e+2jets+X'+Y 



HI preliminary 
NLO, R-1.00 - 
NLO. R=0.62 



10 11 12 
E J T e "[GeV] 



2(10 
180 
160 
140 
E 120 

V-ioo 

1 80 
60 
40 
20 



ep -> e+2jets+X'+Y 
HI preliminary 
NLO,R=l.00 
NLO. R=0.62 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



ep — > e+2jets+X'+Y 
HI preliminary 




Fig. 3: Differential cross sections for diffractive dijet photoproduction as measured by HI with high-i5y e cuts and compared 
to NLO QCD without (R = 1) and with (R = 0.62) global suppression (color online). 
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Fig. 4: Differential cross sections for diffractive dijet photoproduction as measured by HI with high-B^f' cuts and compared 
to NLO QCD with global, resolved, and resolved/direct-IS suppression. 



tionof E J T etl , x° bs and z°p s . For the low-i^f cuts, the resulting suppression factor is R = 0.46 db 0.14, 
which gives in the lowest Ej? -bin a cross section equal to the experimental data point. The error comes 
from the combined experimental statistical and systematic error. The theoretical error due to the scale 
variation is taken into account when comparing to the three distributions. The results of this comparison 
are shown in Figs. la-c. With the exception of Fig. la , where the comparison of da/dE^ 11 is shown, 
the other two plots are such that the data points lie outside the error band based on the scale variation for 
the unsuppressed case. However, the predictions with suppression R = 0.46 agree nicely with the data 
inside the error bands from the scale variation. Most of the data points even agree with the R = 0.46 
predictions inside the much smaller experimental errors. In da/dE^ 11 (see Fig. la) the predictions for 
the second and third bins lie outside the data points with their errors. For R = 1 and R = 0.46 this cross 
sections falls off stronger with increasing Ei^ tl than the data, the normalization being of course about 
two times larger for R = 1. In particular, the third data point agrees with the R = 1 prediction. This 
means that the suppression decreases with increasing E 3 ^ tl (see also Fig. [5] below). This behavior was 
already apparent when we analyzed the first preliminary HI data [12, 13]. Such a behavior points in the 
direction that a suppression of the resolved cross section only would give better agreement with the data, 
as we shall see below. The survival probability R = 0.46 ± 0.14 agrees with the result in [14], which 
quotes R = 0.51 ± 0.01 (stat.) ± 0.10 (syst.), determined by fitting the integrated cross section. From 
our comparison we conclude that the low-^ et data show a global suppression of the order of two in 
complete agreement with the results [12, 13] and [9] based on earlier preliminary and final HI data [9]. 
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Fig. 5: Ratio of the Ei^ a -distribution as measured by HI with low-E^ (left) and high-E 3 T et cuts (right) to the NLO QCD 
prediction without (full), with resolved-only (dashed), and with additional direct IS suppression (dotted). 



Next we want to answer the second question, whether the data could be consistent with a sup- 
pression of the resolved component only. For this purpose we have calculated the cross sections in two 
additional versions: (i) suppression of the resolved cross section and (ii) suppression of the resolved 
cross section plus the NLO direct part which depends on the factorization scale at the photon vertex [4]. 
The suppression factors needed for the two versions will, of course, be different. We determine them 
again by fitting the measured da/dE^ 11 for the lowest £'^f* 1 -bin (see Fig. 2a). Then, the suppression 
factor for version (i) is R = 0.35 (denoted res in the figures), and for version (ii) it is R = 0.32 (de- 
noted res+dir-IS). The results for da/dE 3 T etl , da/dx° bs and da/dzp are shown in Figs. 2a-c, while the 
six other distributions can be found in [15]. We also show the global (direct and resolved) suppression 
prediction with R = 0.46 already shown in Figs. la-c. For the cross section as a function of Zp s , the 
agreement with the global suppression (R = 0.46) and the resolved suppression (R = 0.35 or R = 0.32) 
is comparable. For da / dE^ 11 , the agreement improves considerably for the resolved suppression only 
(note the logarithmic scale in Fig. 2a). The global suppression factor could, of course, be E^-dependent, 
although we see no theoretical reason for such a dependence. For da/dx° bs , which is usually considered 
as the characteristic distribution for distinguishing global versus resolved suppression, the agreement 
with resolved suppression does not improve. Unfortunately, this cross section has the largest hadronic 
corrections of the order of (25 — 30)% [14]. Second, also for the usual photoproduction of dijets the 
comparison between data and theoretical results has similar problems in the large x° bs -bin [19], although 
the .E^-cut is much larger there. In total, we are tempted to conclude from the comparisons in Figs. 
2a-c that the predictions with a resolved-only (or resolved+direct-IS) suppression are consistent with the 
new low-.Ep 6 * HI data [14] and the survival probability is R = 0.35 (only resolved suppression) and 
R = 0.32 (resolved plus direct-IS suppression), respectively. 

The same comparison of the high-E^' data of HI [14] with the various theoretical predictions is 
shown in the following figures. The global suppression factor is obtained again from a fit to the smallest 
£^ e<1 -bin. It is equal to R = 0.62 ± 0.16, again in agreement with the HI result R = 0.62 ± 0.03 (stat.) 
± 0.14 (syst.) [14]. The same cross sections as for the low-E^f* comparison are shown in Figs. 3a-c for 
the two cases R = 1 (no suppression) and R = 0.62 (global suppression), while the six others can again 
be found in [15]. As before with the exception of da/dE 3 ^ 11 and dcr/dM\2 (not shown), most of the 
data points lie outside the R = 1 results with their error bands and agree with the suppressed prediction 
with R = 0.62 inside the respective errors. However, compared to the results in Figs, la-c the distinction 
between the R = 1 band and the R = 0.62 band and the data is somewhat less pronounced, which is 
due to the larger suppression factor. We also tested the prediction for the resolved (resolved+direct-IS) 
suppression, which is shown in Figs. 4a-c. The suppression factor fitted to the smallest bin came out 
as R = 0.38 (res) and R = 0.30 (res+dir-IS), which are almost equal to the corresponding suppression 



ep -> e+2jets+X'+Y 



ZEUS 99-00 x 
NLO, R= 1.00 
NLO, R-0.71 



10 11 12 13 14 15 
Ef'lCeV] 



ep — > e+2jets+X'+Y 

< 7.F.US 99-00 x 1 .15 
NLO. R=l 00 
NT.O, R=071 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



ep -> e+2jets+X'+Y 












— - 


• 








ZEUS 99-00 X 1.15 
NLO, R=1.00 
NLO,R=0.71 





0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Fig. 6: Differential cross sections for diffractive dijet photoproduction as measured by ZEUS and compared to NLO QCD 
without (R — 1) and with (R — 0.71) global suppression (color online). 
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Fig. 7: Differential cross sections for diffractive dijet photoproduction as measured by ZEUS and compared to NLO QCD with 
global, resolved, and resolved/direct-IS suppression. 



factors derived from the low-Ej, data. In most of the comparisons it is hard to observe any preference for 
the global against the pure resolved (resolved plus direct-IS) suppression. We remark that the suppression 
factor for the global suppression is increased by 35%, if we go from the low-E^ to the high-E^ 6 * data, 
whereas for the resolved suppression this increase is only 9%. Under the assumption that the suppression 
factor should not depend on Ej* 1 , we would conclude that the resolved suppression would be preferred, 
as can also be seen from Fig. [5] A global suppression is definitely observed also in the high-E^f* data 
and the version with resolved suppression explains the data almost as well as with the global suppression. 

In Fig. [5] we show the ratio of of the E 3 ^ 1 -distribution as measured by HI to the NLO QCD 
prediction without (full), with resolved-only (dashed), and with additional direct IS suppression (dotted). 
Within the experimental errors, obviously only the former, but not the latter are E^-dependent. 



3 Comparison with ZEUS data 

In this section we shall compare our predictions with the final analysis of the ZEUS data, which was 
published this year [10], in order to see whether they are consistent with the large-£^ e< data of HI. The 
kinematic cuts [10] are almost the same as in the high-EV^f* HI measurements. The only major difference 
to the HI cuts is the larger range in the variable y. Therefore the ZEUS cross sections will be larger than 
the corresponding HI cross sections. The constraint on My is not explicitly given in the ZEUS publica- 
tion [10]. They give the cross section for the case that the diffractive final state consists only of the proton. 
For this they correct their measured cross section by subtracting in all bins the estimated contribution of 



a proton-dissociative background of 16%. When comparing to the theoretical predictions they multiply 
the cross section with the factor 0.87 in order to correct for the proton-dissociative contributions, which 
are contained in the DPDFs 'HI 2006 fit A' and 'HI 2006 fit B' by requiring My < 1.6 GeV. We do 
not follow this procedure. Instead we leave the theoretical cross sections unchanged, i.e. they contain a 
proton-dissociative contribution with My < 1.6 GeV and multiply the ZEUS cross sections by 1.15 to 
include the proton-dissociative contribution. This means that the so multiplied ZEUS cross sections have 
the same proton dissociative contribution as is in the DPDF fits of HI [8]. Since the ZEUS collaboration 
did measurements only for the high-E 1 ^ 6 * cuts, E^ 1 ^ 2 ' > 7.5 (6.5) GeV, we can only compare to those. 
In this comparison we shall follow the same strategy as before. We first compared to the predictions with 
no suppression (R = 1) and then determine a suppression factor by fitting da/dE 3 ^ 11 to the smallest 
E^f -bin. Then we compared to the cross sections as a function of the seven observables x° bs , Zjp s , xjp, 
E^* , y, Mx and rf etl instead of the nine variables in the HI analysis. The distribution in y is equiva- 
lent to the M^-distribution in [14]. The theoretical predictions for these differential cross sections with no 
suppression factor (R = 1) are shown in Figs. 6a-g of [1 1], together with their scale errors and compared 
to the ZEUS data points, and a selection is shown in Fig. [6] Except for the x° bs - and E 3 ^ tl -distributions, 
most of the data points lie outside the theoretical error bands for R = 1. In particular, in Figs. 6b, c, 
e, f and g, most of the points lie outside. This means that most of the data points disagree with the 
unsuppressed prediction. Next, we determine the suppression factor from the measured da / dE 3 ^ tl at the 
lowest E^ tl -bm, 7.5 GeV < E^f* < 9.5 GeV, and obtain R = 0.71. This factor is larger by a factor of 
1.15 than the suppression factor from the analysis of the high-Ejf* data from HI. Curiously, this factor 
is exactly equal to the correction factor we had to apply to restore the dissociative proton contribution. 
Taking the total experimental error of ±7% from the experimental cross section da/dE 3 ^ 11 in the first 
bin into account, the ZEUS suppression factor is 0.71 ± 0.05 to be compared to 0.62 ± 0.14 in the HI 
analysis [14], so that both suppression factors agree inside the experimental errors. 

If we now check how the predictions for R = 0.71 compare to the data points inside the theoretical 
errors, we observe from Figs. 6a-g of Ref . [11] that, with the exception of da / dz°p s and da / dE 3 ^ 11 , most 
of the data points agree with the predictions. This is quite consistent with the HI analysis (see above) 
and leads to the conclusion that also the ZEUS data agree much better with the suppressed predictions 
than with the unsuppressed prediction. In particular, the global suppression factor agrees with the global 
suppression factor obtained from the analysis of the HI data inside the experimental error. 

Similarly as in the previous section we compared the ZEUS data also with the assumption that 
the suppression results only from the resolved cross section. Here, we consider again (i) only resolved 
suppression (res) and (ii) resolved plus direct suppression of the initial-state singular part (res+dir-IS). 
For these two models we obtain the suppression factors R = 0.53 and R = 0.45, respectively, where 
these suppression factors are again obtained by fitting the data point at the first bin of da / dE 3 ^ tl . The 
comparison to the global suppression with R = 0.71 and to the data is shown in Figs. 7a-g of [11] and 
a selection in Fig. [7] In general, we observe that the difference between global suppression and resolved 
suppression is small, i.e. the data points agree with the resolved suppression as well as with the global 
suppression. 

4 Conclusion 

In conclusion, we found that most of the data points of diffractive dijet photoproduction in the latest HI 
analyses with low- and high-£^ e< cuts and in the final ZEUS analysis with the same high-E'jf* cuts dis- 
agree with NLO QCD predictions within experimental and theoretical errors. When global factorization 
breaking is assumed in both the direct and resolved contributions, the resulting suppression factor would 
have to be E^ -dependent, although we see no theoretical motivation for this assumption. Suppressing 
only the resolved or in addition the direct initial-state singular contribution by about a factor of three, as 
motivated by the proof of factorization in point-like photon-hadron scattering and predicted by absorp- 



tive models [7], the agreement between theory and data is at least as good as for global suppression, and 
no Ej^ 1 -dependence of the survival probability is needed. 
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